Intercomparison of Nine Micrometeorological Stations during the BEAREX08 Field Campaign

Land–atmosphere interactions play a critical role in regulating numerous meteorological, hydrological, and environmental processes. Investigating these processes often requires multiple measurement sites representing a range of surface conditions. Before these measurements can be compared, however, it is imperative that the differences among the instrumentation systems are fully characterized. Using data collected as a part of the 2008 Bushland Evapotranspiration and Agricultural Remote Sensing Experiment (BEAREX08), measurements from nine collocated eddy covariance (EC) systems were compared with the twofold objective of 1) characterizing the interinstrument variation in the measurements, and 2) quantifying the measurement uncertainty associated with each system. Focusing on the three turbulent fluxes (heat, water vapor, and carbon dioxide), this study evaluated the measurement uncertainty using multiple techniques. The results of the analyses indicated that there could be substantial variability in the uncertainty estimates because of the advective conditions that characterized the study site during the afternoon and evening hours. However, when the analysis was limited to nonadvective, quasi-normal conditions, the response of the nine EC stations were remarkably similar. For the daytime period, both the method of Hollinger and Richardson and the method of Mann and Lenschow indicated that the uncertainty in the measurements of sensible heat, latent heat, and carbon dioxide flux were approximately 13 W m‒2, 27 W m‒2, and 0.10 mg m‒2 s‒1, respectively. Based on the results of this study, it is clear that advection can greatly increase the uncertainty associated with EC flux measurements. Since these conditions, as well as other phenomena that could impact the measurement uncertainty, are often intermittent, it may be beneficial to conduct uncertainty analyses on an ongoing basis

Influence of wind direction on the surface roughness of vineyards

Remote sensing-based models are the most viable means of collecting the high-resolution spatially distributed estimates of evaporative water loss needed to manage irrigation and ensure the effective use of limited water resources. However, due to the unique canopy structure and configuration of vineyards, these models may not be able to adequately describe the physical processes driving evapotranspiration from vineyards. Using data collected from 2014 to 2016 as a part of the Grape Remote sensing Atmospheric Profile and Evapotranspiration Experiment (GRAPEX), the twofold objective of this study was to (1) identify the relationship between the roughness parameters, zero-plane displacement height (do) and roughness length for momentum (zo), and local environmental conditions, specifically wind direction and vegetation density and (2) determine the effect of using these relationships on the ability of the remote sensing-based Two-Source Energy Balance (TSEB) model to estimate the sensible (H) and latent (λE) heat fluxes. Although little variation in do was identified during the growing season, a well-defined sigmoidal relationship was observed between zo and wind direction. When the output from a version of the TSEB model incorporating these relationships (TSEBVIN) was compared to output from the standard model (TSEBSTD), there were large changes to the roughness parameters, particularly zo, but only modest changes in the turbulent fluxes. When the output from TSEBVIN was compared to that of a version using a parameterization scheme representing open canopies (TSEBOPN), the mean absolute difference between the estimates of do and zo were 0.44 m and 0.25 m, respectively. While these values represent differences in excess of 45%, the turbulent fluxes differed by just 13 W m−2 or 10%, on average. The results suggest that the TSEB model is largely insensitive to changes in the roughness parameters for the range in roughness values evaluated in this study. This also suggests that the requirement for highly accurate roughness values has limited utility in the application of the TSEB model in vineyard systems. Since there is no significant advantage to using the more complex TSEBOPN and TSEBVIN models, it is recommended that the standard model be used.info:eu-repo/semantics/acceptedVersio

Field-scale mapping of evaporative stress indicators of crop yield: An application over Mead, NE, USA

The Evaporative Stress Index (ESI) quantifies temporal anomalies in a normalized evapotranspiration (ET) metric describing the ratio of actual-to-reference ET (fRET) as derived from satellite remote sensing. At regional scales (3–10 km pixel resolution), the ESI has demonstrated the capacity to capture developing crop stress and impacts on regional yield variability in water-limited agricultural regions. However, its performance in some regions where the vegetation cycle is intensively managed appears to be degraded due to spatial and temporal limitations in the standard ESI products. In this study, we investigated potential improvements to ESI by generating maps of ET, fRET, and fRET anomalies at high spatiotemporal resolution (30-m pixels, daily time steps) using a multi-sensor data fusion method, enabling separation of landcover types with different phenologies and resilience to drought. The study was conducted for the period 2010–2014 covering a region around Mead, Nebraska that includes both rainfed and irrigated crops. Correlations between ESI and measurements of maize yield were investigated at both the field and county level to assess the potential of ESI as a yield forecasting tool. To examine the role of crop phenology in yield-ESI correlations, annual input fRET time series were aligned by both calendar day and by biophysically relevant dates (e.g. days since planting or emergence). At the resolution of the operational U.S. ESI product (4 km), adjusting fRET alignment to a regionally reported emergence date prior to anomaly computation improves r2 correlations with county-level yield estimates from 0.28 to 0.80. At 30-m resolution, where pure maize pixels can be isolated from other crops and landcover types, county-level yield correlations improved from 0.47 to 0.93 when aligning fRET by emergence date rather than calendar date. Peak correlations occurred 68 days after emergence, corresponding to the silking stage for maize when grain development is particularly sensitive to soil moisture deficiencies. The results of this study demonstrate the utility of remotely sensed ET in conveying spatially and temporally explicit water stress information to yield prediction and crop simulation models

Mapping daily evapotranspiration at Landsat spatial scales during the BEAREX’08 field campaign

Robust spatial information about environmental water use at field scales and daily to seasonal timesteps will benefit many applications in agriculture and water resource management. This information is particularly critical in arid climates where freshwater resources are limited or expensive, and groundwater supplies are being depleted at unsustainable rates to support irrigated agriculture as well as municipal and industrial uses. Gridded evapotranspiration (ET) information at field scales can be obtained periodically using land–surface temperature-based surface energy balance algorithms applied to moderate resolution satellite data from systems like Landsat, which collects thermal-band imagery every 16 days at a resolution of approximately 100 m. The challenge is in finding methods for interpolating between ET snapshots developed at the time of a clear-sky Landsat overpass to provide complete daily time-series over a growing season. This study examines the efficacy of a simple gap-filling algorithm designed for applications in data-sparse regions, which does not require local ground measurements of weather or rainfall, or estimates of soil texture. The algorithm relies on general conservation of the ratio between actual ET and a reference ET, generated from satellite insolation data and standard meteorological fields from a mesoscale model. The algorithm was tested with ET retrievals from the Atmosphere–Land Exchange Inverse (ALEXI) surface energy balance model and associated DisALEXI flux disaggregation technique, which uses Landsat-scale thermal imagery to reduce regional ALEXI maps to a finer spatial resolution. Daily ET at the Landsat scale was compared with lysimeter and eddy covariance flux measurements collected during the Bushland Evapotranspiration and Agricultural Remote sensing EXperiment of 2008 (BEAREX08), conducted in an irrigated agricultural area in the Texas Panhandle under highly advective conditions. The simple gap-filling algorithm performed reasonably at most sites, reproducing observed cumulative ET to within 5–10% over the growing period from emergence to peak biomass in both rainfed and irrigated fields

Evaporative loss from irrigated interrows in a highly advective semi-arid agricultural area

Agricultural productivity has increased in the Texas High Plains at the cost of declining water tables, putting at risk the sustainability of the Ogallala Aquifer as a principal source of water for irrigated agriculture. This has led area producers to seek alternative practices that can increase water use efficiency (WUE) through more careful management of water. One potential way of improving WUE is by reducing soil evaporation (E), thus reducing overall evapotranspiration (ET). Before searching for ways to reduce E, it is first important to quantify E and understand the factors that determine its magnitude. The objectives of this study were (1) to quantify E throughout part of the growing season for irrigated cotton in a strongly advective semi-arid region; (2) to study the effects of LAI, days after irrigation, and measurement location within the row on the E/ET fraction; and (3) to study the ability of microlysimeter (ML) measures of E combined with sap flow gage measures of transpiration (T) to accurately estimate ET when compared with weighing lysimeter ET data and to assess the E/T ratio. The research was conducted in an irrigated cotton field at the Conservation & Production Research Laboratory of the USDA-ARS, Bushland, TX. ET was measured by a large weighing lysimeter, and E was measured by 10 microlysimeters that were deployed in two sets of 5 across the interrow. In addition, 10 heat balance sap flow gages were used to determine T. A moderately good agreement was found between the sum E + T and ET (SE = 1 mm or ~10% of ET). It was found that E may account for \u3e50% of ET during early stages of the growing season (LAI \u3c 0.2), significantly decreasing with increase in LAI to values near 20% at peak LAI of three. Measurement location within the north-south interrows had a distinct effect on the diurnal pattern of E, with a shift in time of peak E from west to east, a pattern that was governed by the solar radiation reaching the soil surface. However, total daily E was unaffected by position in the interrow. Under wet soil conditions, wind speed and direction affected soil evaporation. Row orientation interacted with wind direction in this study such that aerodynamic resistance to E usually increased when wind direction was perpendicular to row direction; but this interaction needs further study because it appeared to be lessened under higher wind speeds

Mapping evapotranspiration with high-resolution aircraft imagery over vineyards using one- and two-source modeling schemes

Thermal and multispectral remote sensing data from low-altitude aircraft can provide high spatial resolution necessary for sub-field ( 10 m) and plant canopy (1 m) scale evapotranspiration (ET) monitoring. In this study, highresolution (sub-meter-scale) thermal infrared and multispectral shortwave data from aircraft are used to map ET over vineyards in central California with the two-source energy balance (TSEB) model and with a simple model having operational immediate capabilities called DATTUTDUT (Deriving Atmosphere Turbulent Transport Useful To Dummies Using Temperature). The latter uses contextual information within the image to scale between radiometric land surface temperature (TR) values representing hydrologic limits of potential ET and a non-evaporative surface. Imagery from 5 days throughout the growing season is used for mapping ET at the sub-field scale. The performance of the two models is evaluated using tower-based measurements of sensible (H) and latent heat (LE) flux or ET. The comparison indicates that TSEB was able to derive reasonable ET estimates under varying conditions, likely due to the physically based treatment of the energy and the surface temperature partitioning between the soil/cover crop inter-row and vine canopy elements. On the other hand, DATTUTDUT performance was somewhat degraded presumably because the simple scaling scheme does not consider differences in the two sources (vine and inter-row) of heat and temperature contributions or the effect of surface roughness on the efficiency of heat exchange. Maps of the evaporative fraction (EFDLE/(H CLE)) from the two models had similar spatial patterns but different magnitudes in some areas within the fields on certain days. Large EF discrepancies between the models were found on 2 of the 5 days (DOY 162 and 219) when there were significant differences with the tower-based ET measurements, particularly using the DATTUTDUT model. These differences in EF between the models translate to significant variations in daily water use estimates for these 2 days for the vineyards. Model sensitivity analysis demonstrated the high degree of sensitivity of the TSEB model to the accuracy of the TR data, while the DATTUTDUT model was insensitive to systematic errors in TR as is the case with contextual-based models. However, it is shown that the study domain and spatial resolution will significantly influence the ET estimation from the DATTUTDUT model. Future work is planned for developing a hybrid approach that leverages the strengths of both modeling schemes and is simple enough to be used operationally with high-resolution imagery

To What Extent Does the Eddy Covariance Footprint Cutoff Influence the Estimation of Surface Energy Fluxes Using Two Source Energy Balance Model and High-Resolution Imagery in Commercial Vineyards?

Validation of surface energy fluxes from remote sensing sources is performed using instantaneous field measurements obtained from eddy covariance (EC) instrumentation. An eddy covariance measurement is characterized by a footprint function / weighted area function that describes the mathematical relationship between the spatial distribution of surface flux sources and their corresponding magnitude. The orientation and size of each flux footprint / source area depends on the micro-meteorological conditions at the site as measured by the EC towers, including turbulence fluxes, friction velocity (ustar), and wind speed, all of which influence the dimensions and orientation of the footprint. The total statistical weight of the footprint is equal to unity. However, due to the large size of the source area / footprint, a statistical weight cutoff of less than one is considered, ranging between 0.85 and 0.95, to ensure that the footprint model is located inside the study area. This results in a degree of uncertainty when comparing the modeled fluxes from remote sensing energy models (i.e., TSEB2T) against the EC field measurements. In this research effort, the sensitivity of instantaneous and daily surface energy flux estimates to footprint weight cutoffs are evaluated using energy balance fluxes estimated with multispectral imagery acquired by AggieAir sUAS (small Unmanned Aerial Vehicle) over commercial vineyards near Lodi, California, as part of the ARS-USDA Agricultural Research Service’s Grape Remote Sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX) project. The instantaneous fluxes from the eddy covariance tower will be compared against instantaneous fluxes obtained from different TSEB2T aggregated footprint weights (cutoffs). The results indicate that the size, shape, and weight of pixels inside the footprint source area are strongly influenced by the cutoff values. Small cutoff values, such as 0.3 and 0.35, yielded high weights for pixels located within the footprint domain, while large cutoffs, such as 0.9 and 0.95, result in low weights. The results also indicate that the distribution of modelled LE values within the footprint source area are influenced by the cutoff values. A wide variation in LE was observed at high cutoffs, such as 0.90 and 0.95, while a low variation was observed at small cutoff values, such as 0.3. This happens due to the large number of pixel units involved inside the footprint domain when using high cutoff values, whereas a limited number of pixels are obtained at lower cutoff values

Estimation of Evapotranspiration and Energy Fluxes Using a Deep-Learning-Based High-Resolution Emissivity Model and the Two-Source Energy Balance Model with sUAS Information

Surface temperature is necessary for the estimation of energy fluxes and evapotranspiration from satellites and airborne data sources. For example, the Two-Source Energy Balance (TSEB) model uses thermal information to quantify canopy and soil temperatures as well as their respective energy balance components. While surface (also called kinematic) temperature is desirable for energy balance analysis, obtaining this temperature is not straightforward due to a lack of spatially estimated narrowband (sensor-specific) and broadband emissivities of vegetation and soil, further complicated by spectral characteristics of the UAV thermal camera. This study presents an effort to spatially model narrowband and broadband emissivities for a microbolometer thermal camera at UAV information resolution (~0.15 m) based on Landsat and NASA HyTES information using a deep learning (DL) model. The DL model is calibrated using equivalent optical Landsat / UAV spectral information to spatially estimate narrowband emissivity values of vegetation and soil in the 7–14- nm range at UAV resolution. The resulting DL narrowband emissivity values were then used to estimate broadband emissivity based on a developed narrowband-broadband emissivity relationship using the MODIS UCSB Emissivity Library database. The narrowband and broadband emissivities were incorporated into the TSEB model to determine their impact on the estimation of instantaneous energy balance components against ground measurements. The proposed effort was applied to information collected by the Utah State University AggieAir small Unmanned Aerial Systems (sUAS) Program as part of the ARS-USDA GRAPEX Project (Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment) over a vineyard located in Lodi, California. A comparison of resulting energy balance component estimates, with and without the inclusion of high-resolution narrowband and broadband emissivities, against eddy covariance (EC) measurements under different scenarios are presented and discussed

Implications of Soil and Canopy Temperature Uncertainty in the Estimation of Surface Energy Fluxes Using TSEB2T and High-Resolution Imagery in Commercial Vineyards

Estimation of surface energy fluxes using thermal remote sensing–based energy balance models (e.g., TSEB2T) involves the use of local micrometeorological input data of air temperature, wind speed, and incoming solar radiation, as well as vegetation cover and accurate land surface temperature (LST). The physically based Two-source Energy Balance with a Dual Temperature (TSEB2T) model separates soil and canopy temperature (Ts and Tc) to estimate surface energy fluxes including Rn, H, LE, and G. The estimation of Ts and Tc components for the TSEB2T model relies on the linear relationship between the composite land surface temperature and a vegetation index, namely NDVI. While canopy and soil temperatures are controlling variables in the TSEB2T model, they are influenced by the NDVI threshold values, where the uncertainties in their estimation can degrade the accuracy of surface energy flux estimation. Therefore, in this research effort, the effect of uncertainty in Ts and Tc estimation on surface energy fluxes will be examined by applying a Monte Carlo simulation on NDVI thresholds used to define canopy and soil temperatures. The spatial information used is available from multispectral imagery acquired by the AggieAir sUAS Program at Utah State University over vineyards near Lodi, California as part of the ARS-USDA Agricultural Research Service’s Grape Remote Sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX) project. The results indicate that LE is slightly sensitive to the uncertainty of NDVIs and NDVIc. The observed relative error of LE corresponding to NDVIs uncertainty was between -1% and 2%, while for NDVIc uncertainty, the relative error was between -2.2% and 1.2%. However, when the combined NDVIs and NDVIc uncertainties were used simultaneously, the domain of the observed relative error corresponding to the absolute values of |ΔLE| was between 0% and 4%

Incorporation of Unmanned Aerial Vehicle (UAV) Point Cloud Products into Remote Sensing Evapotranspiration Models

In recent years, the deployment of satellites and unmanned aerial vehicles (UAVs) has led to production of enormous amounts of data and to novel data processing and analysis techniques for monitoring crop conditions. One overlooked data source amid these efforts, however, is incorporation of 3D information derived from multi-spectral imagery and photogrammetry algorithms into crop monitoring algorithms. Few studies and algorithms have taken advantage of 3D UAV information in monitoring and assessment of plant conditions. In this study, different aspects of UAV point cloud information for enhancing remote sensing evapotranspiration (ET) models, particularly the Two-Source Energy Balance Model (TSEB), over a commercial vineyard located in California are presented. Toward this end, an innovative algorithm called Vegetation Structural-Spectral Information eXtraction Algorithm (VSSIXA) has been developed. This algorithm is able to accurately estimate height, volume, surface area, and projected surface area of the plant canopy solely based on point cloud information. In addition to biomass information, it can add multi-spectral UAV information to point clouds and provide spectral-structural canopy properties. The biomass information is used to assess its relationship with in situ Leaf Area Index (LAI), which is a crucial input for ET models. In addition, instead of using nominal field values of plant parameters, spatial information of fractional cover, canopy height, and canopy width are input to the TSEB model. Therefore, the two main objectives for incorporating point cloud information into remote sensing ET models for this study are to (1) evaluate the possible improvement in the estimation of LAI and biomass parameters from point cloud information in order to create robust LAI maps at the model resolution and (2) assess the sensitivity of the TSEB model to using average/nominal values versus spatially-distributed canopy fractional cover, height, and width information derived from point cloud data. The proposed algorithm is tested on imagery from the Utah State University AggieAir sUAS Program as part of the ARS-USDA GRAPEX Project (Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment) collected since 2014 over multiple vineyards located in California. The results indicate a robust relationship between in situ LAI measurements and estimated biomass parameters from the point cloud data, and improvement in the agreement between TSEB model output of ET with tower measurements when employing LAI and spatially-distributed canopy structure parameters derived from the point cloud data